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FIELD OF THE INVENTION 

The present invention relates to a method and to experimental devices for 
evaluating, with a single equipment and simvdtaneously, the porosity and the absolute 
permeability of any fragmented natural or artificial more or less porous medium from 
5 rock samples taken from this medium. The samples can be for example rocks from a 
zone of an underground fluid reservoir (hydrocarbon reservoir, aquifer, etc.) or possibly 
rocks of very low permeability such as clays or shales. The fragments are for example 
obtained during well drilling operations : drill cuttings or fragments obtained by 
crushing bigger samples : cores or lateral cores taken in a well. 

10 The current petroleum context leads operators to take an interest in new zones (deep 

offshore) and in new reservoir types (marginal structures located close to existing 
surface installations). Considering the drilling costs linked with the difiHcult 
environment of these new discoveries or with the limited size of certain stmctures, 
operators can no longer afford to drill complementary appraisal wells without risking 

15 compromising the economic viability of the project. The development strategy set 
before starting production is therefore less strict so as to allow « real-time » adaptation 
to the nature of the information collected by drilling production wells. This is referred to 
as appraisal development. 

Petrophysical measurements have a key role in the assessment of the quality of a 
20 reservoir. However, the time required for this type of measurements is often very long 
and therefore incompatible with the reactivity necessary for the success of appraisal 
developments. New, faster and less expensive evaluation means are thus required as a 
support to decision-making. 
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Cuttings carried along by the drillinig mud have been subjected to examinations in 
the field for a long time. They are carried out by the crews in charge of mud logging and 
they essentially help to complete the description of the geologic layers drilled, 
performed from logs. 

5 BACKGROUND OF THE INVENTION 

Work has already been done to try to evaluate petrophysical properties from drill 
cuttings. The acoustic properties as regards S and P waves (shear and compressional 
waves) have been measured for example. Various parameters have also been studied, 
such as the hardness and the deformation of the rock fragments, or their porosity and 
10 permeability. 

According to a first known method intended for permeability measurement, the rock 
fragment is first coated with resin. A thin slice is cut out of the coated rock and placed 
in a measuring cell. It comprises means for injecting a fluid under pressure at a 
controlled flow rate and means for measuring the pressure drop created by the sample. 
15 Since the resin is impermeable, the absolute permeability is deduced from Darcy's 
equation by taking account of the real surface area occupied by the rock fragments. 

This method is for example described by : 

- Santarelli F.J., et al. ; « Formation evaluation from logging on cuttings », SPERE, June 
1998, or 

20 - Marsala A.F., et al. ; « Transient Method Implemented under Unsteady State 
Conditions for Low and Very Low Permeability Measurements on Cuttings », 
SPE/ISRM No.47202, Trondheim, 8-10 July 1998. 
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This type of measurement can only be obtained in the laboratory after long 
fragment conditioning operations. 

Another well-known method is based on a NMR (Nuclear Magnetic Resonance) 
measurement which is directly carried out on the cuttings after prior washing followed 
5 by brine saturation. This type of measurement gives a porosity value that can be directly 
used. The permeability K is determined by means of correlations of the same nature as 
those used within the context of NMR logging. 

An illustration of this method can be found in the following document : 

- Nigh E., et al. ; P-K™ : Wellsite Determination of Porosity and Permeability Using 
10 Drilling Cuttings, GWLS Journal, Vol, 13, No. 1, Dec. 1984. 

Patent application EP- 1,167,948 arid patent applications FR-02/02,242 and 
03/00,429 describe several systems for evaluating physical parameters such as the 
absolute permeability of porous rocks from a zone of an underground reservoir, from 
rock fragments taken in this zone, such as rock fragments carried along by the drilling 
15 mud. It comprises a stage of immersing fragments contained in a containment chamber 
in a viscous fluid and of placing the chamber containing the Augments under pressure 
according to several implementation modes so as to compress the gas trapped in the 
pores of the rock. 

According to an implementation mode, a fluid is injected at a pressure that increases 
20 with time, up to a predetennined pressure threshold, so as to compress the gas trapped 
in the pores of the rock. This injection stage is followed by a relaxation stage where 
injection is stopped. The pressure variation during these two successive stages is 
recorded. 
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According to another implementation mode, the chamber containing the fragments 
is communicated with a vessel containing a viscous fluid at a predetermined pressure. 
The communication stage being short, a &st pressure increase in the chamber and 
compression of the gas trapped in the pores of the rock are obtained, followed by a 
5 relaxation period after isolation of the chamber, and the pressure evolution in the 
chamber during the two periods is measured. 

According to another implementation mode, the chamber containing the fragments 
is conmiunicated with a vessel containing a viscous fluid at a predetermined constant 
pressure so as to cause a fast and prolonged pressure increase in the chamber and 
10 compression of the gas trapped in the pores of the rock, and the evolution of the volume 
of fluid injected is measured as a function of time. 

Whatever the operating mode, the evolution of the pressure or of the volume of 
fluid injected having been modelled from initial values selected for the physical 
parameters of the fragments, the computer adjusts them iteratively so as to obtain the 
15 best possible match between the modelled pressure or volume curve and the pressure or 
volume curve really measured. 

SUMMARY OF THE INVENTION 

The object of the method according to the invention is to evaluate, with a single 
equipment, physical parameters such as the absolute permeability and the porosity of a 
20 fragmented natural or artificial porous medium from rock fragments taken in this 
medium. It comprises a stage of immersing fragments contained in a containment 
chamber in a (more or less viscous) liquid and of intermittently communicating the 
chamber containing the fragments with a source of liquid under pressure comprising an 
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accumulator containing an elastic volume of liquid so as to compress the gas trapped in 
the pores of the rock, a stage of measuring the volume of liquid injected, a stage of 
modelling the evolution of the volume of liquid injected from a priori selected initial 
values at least for the permeability and the residual gas saturation, and a stage of 
5 iterative adjustment of the value of the physical parameters of the rock fragments so as 
to obtain the best possible adjustment of the modelled evolution of the volume injected 
with the measured evolution of the volume injected in the chamber. 

The method is characterized in that, during the stage of immersing in the liquid the 
fragments contained in the containment chamber, the volume of liquid injected is 
10 measured by measuring the concomitant variation of the pressure prevailing in said 
accumulator. 

The modelling stage can also be carried out from a priori selected initial porosity 
values. It can also comprise a prior stage of feeding the washed and dried rock 
fragments into the containment chamber which is first communicated with a gas tank at 
15 a predetermined pressure, so . as to determine the solid volume of said fragments, with 
measurement of the envelope volume and of the mass of the fragments so as to deduce 
the porosity and the density of the rock fragments. In this case, the modelling stage 
comprises modelling the evolution of the volume of liquid injected from a priori 
. selected initial permeability and residual gas saturation values, and from the value 
20 previously measured for the porosity. 

According to an implementation mode, the stage of communicating the chamber 
with the accumulator is carried out so as to cause a fast pressure increase in the chamber 
and compression of the gas trapped in the pores of the rock, followed by a relaxation 
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period after isolation of the chamber, and it comprises measuring the evolution of the 
pressure in the chamber during these two successive stages. 

The containment chamber can for example be filled with drill cuttings, rock 
fi-agments obtained by crushing cores taken in a well, notably cores obtained by 
5 sidewall coring of a well, rock fi-agments invaded by drilling fluids, previously cleaned 
rock fi-agments, etc. 

The device according to the invention allows to evaluate physical parameters such 
as the absolute permeability and the porosity of a fragmented artificial or natural porous 
medium such as a zone of an underground reservoir, fi-om fi-agments taken fi-om this 

10 medium, comprising a processing system, a containment chamber for the fi-agments, an 
injection assembly for injecting a liquid into the chamber so as to fill the chamber 
containing the rock fi-agments, and for carrying out a cycle comprising a stage of liquid 
injection into the chamber, this assembly comprising an accumulator containing an 
elastic volume of liquid and means controlled by the processing system for controlling 

15 communication of the accumulator with ,the chamber containing the rock fi-agments. 
Means are used to measure the pressure in the chamber. The processing system is suited 
for modelling the evolution of the volume of liquid injected from initial values selected 
for the physical parameters of the rock fragments, and for adjusting iteratively the 
values to be given to these physical parameters so as to obtain the best possible 

20 adjustment between the modelled evolution of the physical quantity and the measured 
evolution of said quantity in the chamber. 

The device is characterized in that it comprises means for measuring pressure 
variations in the accumulator and in that the processing system is suited for computing 
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the volume of liquid injected in the chamber from the accumulator, from the pressure 
variation measured by the pressure measuring means. 

The device comprises for example a gas tank that can be communicated with the 
chamber by means of a valve, a means for measuring the volume envelope (such as a 
5 powder pycnometer for example) for determination of the porosity of the fragments and 
a means for measuring the mass of the fragments. 

The elastic volume of liquid is for example delimited in the accumulator by either a 
gas volume, in which case the relation is obtained by applying the gas law, or by an 
elastic membrane or any other moving element, in which case the relation is obtained by 
10 prior calibration of the deformations or displacements. 

The method proves to be satisfactory for very varied rocks within a wide 
permeability and porosity range. It is very accurate because the measurement of the 
volume of liquid injected is not obtained directly as in the aforementioned patent 
applications but indirectly by measuring the pressure variations in the accumulator, 
15 which can be readily done with high precision. 

Considering the limited space occupied by the device, its ease of use and the 
rapidity with which the measurements and the adjustment between the theoretical data 
and the experimental data can be performed, the method is particularly well-suited for 
field conditions. It is therefore possible to carry out measurements and interpretations 
20 directly in the field in a very short time, which consequently has nothing in conmion 
with the time required to obtain equivalent results by means of laboratory methods. This 
opens up interesting prospects for petrophysical characterization in the petroleum and 
hydrologic sphere : permeability evolution, porosity, rock density, etc., since it is 



possible to take advantage of this new source of information as a support for the 
interpretation of electric logs and to fine down the evaluation of a well in terms of 
production potential. This also opens up important prospects for the petrophysical 
characterization of any other fi^agmented natural or artificial porous medium. 

5 BRIEF DESCRIPTION OF THE FIGURES 

Other features and advantages of the method and of the device according to the 
invention will be clear fi-om reading the description hereafter of embodiments given by 
way of non limitative example, with reference to the accompanying drawings wherein : 

- Figure 1 diagranmiatically shows the device, 

10 - Figure 2 diagrammatically shows the structure of a fi-agment or particle of a porous 
rock wherein the effects of the injection of a high- viscosity liquid such as oil are 
modelled, 

- Figures 3 A and 3B show a pressure evolution measured during a test with the 
corresponding evolution of the volume of liquid injected as a function of time calculated 

15 from the gas law, 

- Figures 4 A to 4C diagrammatically show the curves of the volume of liquid injected as 
a fiinction of time obtained from the pressure evolution in the chamber of the device of 
Fig. 1 for three different rocks, 

- Figure 5 shows the good match between the permeabilities obtained for several rock 
20 types by means of a conventional core measuring method and by means of the method 

according to the invention. 
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- Figures 6A to 6C show, for the previous three rocks, the accuracy obtained in the 
adjustment of the modelled injected volume ounces in relation to the experimental 
curves, 

- Figure 7 shows the good match obtained between the values of the porosity (O or Phic) 
5 of the rock fragments tested and those (Phis) obtained for samples or cores in cases 

where the a priori porosity and permeability values are iteratively adjusted, and 

- Figure 8 shows the result of the porosity determination of the rock fragments tested 
when the porosity (Phic) of the fragments is first measured experimentally. 

DETAILED DESCRIPTION 

10 The device diagrammatically shown in Figure 1 comprises a containment chamber 

1 in which the drill cuttings are initially placed. An injection system 2 conmiunicates 
through a line 3 with the base of a surge tank 4 containing a more or less viscous liquid 
suited to the permeability to be measured, such as oil. The opposite end of surge tank 4 
communicates through a valve V6 with a line LI. A first end of containment chamber 1 

15 communicates with line LI through the agency of two valves VI, V2, The opposite end 
of containment chamber 1 communicates, by means of an isolating valve VS, with a 
separator 6. A manometer 7 is connected to the outlet of chamber 1. The pressure 
variations measured by manometer 1 are acquired by a processor 8 such as a 
microcomputer. Valve V2 is directly controlled by processor 8. Line LI also 

20 communicates by means of a valve V7 with an accumulator 9 containing a more or less 
viscous liquid, placed under a predetermined pressure by a gas cap under pressure, a 
flexible membrane or a moveable element. A manometer 10 is connected to the outlet of 
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accumulator 9. The pressure variations measured thereby are also acquired by processor 
8. 

Fragments available after cleaning, from which all the fluids have been previously 
removed, can also be used. In cases where containment chamber 1 is filled with cleaned 
5 fragments, a gas such as helium for example, coming from a bottle 5, can be injected by 
opening a valve V4 so as to expel the air from the chamber. 

Chamber 1 is also connected, by means of valves VI and V3, to a tank 1 1 filled 
with helium and whose volume is known. Chamber 1 can be communicated with tank 
1 1 initially at a known pressure by opening Vl or V3, valves V2 and V4 being closed. 

10 The device fiirther comprises a balance 12 and a powder pycnometer type 

instrument 13 allowing to measure the envelope volume of the fragments introduced. 

I) Porosity measurement 

Determination of the porosity comprises a stage of acquisition of experimental 
measurements of the envelope volume Ve of the fragmeiits introduced, of the solid 
15 volume Vs of rock introduced and of the mass nu of rock introduced, and a stage of 
calculation of the porosity and of the density of the rock. 

a) Measurement acquisition 

The cleaned and dried fragments are previously weighed on balance 12 and their 
envelope volume is measured by means of instrument 13. The fragments are then fed 
20 into containment chamber 1 placed under gaseous atmosphere (helium for example) 
through conmiunication with tank 5 so as to expel the air. Chamber 1 is then connected 
to tank 1 1 filled with helium and of known volume by opening valves V 1 and V3, 



11 



valves V2 and V4 being closed. The balance pressure allows to deduce the value of the 
solid volume of the rock from the initial pressures in chamber 1 and tank 1 1 and from 
their volumes. All these measurements allow to determine the porosity of the samples. 
The density of the rock is also obtained by measuring the mass of the fragments 
5 introduced. 

Envelope volume Ve is obtained by means of a powder pycnometer according to a 
technique known to specialists. 

Let PI be the initial pressure in 1, Ph the initial pressure in tank 1 1, Pe the balance 
pressure after conmiunication, Ul the volume of chamber 1, Vh the volume of tank 11 
10 and mc the mass of the fragments used. 

Vs, <|) and d are respectively the solid volume, the porosity and the density of the 
rock. We have : 



(Pe-PV) 

^ Ve-Vs 

0 = — 

Ve 

d rock = 

Ve 



Figure 1 1 shows that a very good estimation of the porosity of the rocks tested is 
15 obtained. 



n) Permeability estimation 

Estimation of the absolute permeability essentially comprises three stages : 



1) a stage of acquisition of experimental measurements of the pressure variations 
(or of the evolution of the volume injected) from the drilling fragments, which leads to 
experimental curves, 

2) a stage of modelling the physical phenomena occurring in the drilling fragments 
5 during the same operating cycle, for arbitrary values of physical parameters sought 

(permeability K and porosity <|)) involved in the model, allowing to establish similar 
theoretical curves, and 

3) an adjustment or calibration stage wherein the values to be given to the physical 
parameters involved in the model are determined so as to obtain the best possible 

10 adjustment between the experimental curves and the theoretical curves. 

1) Measurement acquisition 

Chamber 1 is filled with a more or less viscous liquid (with a viscosity suited to the 
permeability to be measured : oil for example, glycerol, water, any aqueous, viscous 
solution, etc.) or by means of the injection system. The liquid occupies the free space 
15 between the drilling fragments and it also flows, by spontaneous imbibition, into the 
rock. A liberation of gas whose intensity and duration depend on the nature of the rock 
(mainly the porosity) takes place. This gas liberation only involves part of the gas. A 
certain residual volume remains trapped in the drilling fragments in form of 
disconnected pockets. 

20 The procedure essentially consists in communicating chamber 1, which contains 

rock fragments C, with accumulator 9 which contains liquid under pressure (Pini) by 
opening valve V2 controlled by control processor 8. 
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Initially, the valve is closed. On the side of accumulator 9, the pressure is Pini, 
whereas on the side of chamber 1, the pressure is equal to the ambient pressure. Valve 
V2 controlled by processor 8 is then opened for some tenths of a second so as to rapidly 
increase the pressure in chamber 1 up to Pmax. The elasticity of the volume of 

5 accumulator 9, obtained by means of a gas cap, of a flexible membrane that deforms 
under the effect of a pressure change or of a moveable element, is previously calibrated 
so as to measure a pressure decrease of some bars in accumulator 9 during the test. This 
pressure change then allows to readily calculate the evolution of the volume of liquid 
injected from accumulator 9 into chamber 1 from the gas law (if a gas cap is used) or 

10 from a calibration curve (if a flexible membrane or a moveable element is used). 

In relation to the procedures described in the aforementioned patent applications, 
this approach allows to simplify the course of the experiment (no flowmeter or 
differential detector to be inserted as previously on line LI connecting the accumulator 
to the chamber to measure the volume injected) while keeping good control of the 
15 volume injected, which facilitates acquisition of the measurements and interpretation of 
the results with the simulator. Besides, elimination of a flowmeter or of a differential 
detector for measuring the volume injected allows : 

- to greatly miniaturize the dimensions of the measuring instrument, so that a 
measurement can be performed with only 1 cm^ of fi-agment, 

20 - to shorten the connecting lines between accumulator 9 and chamber 1, which shortens 
the time during which the chamber containing the fragments is placed under pressure 
and improves the sensitivity of the instrument. It is thus possible to test rocks having an 
average diameter of the order of 1 mm for permeabilities up to several hundred 
milliDarcy, 
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- to measure very small variations of the volume injected (less than 0.03 cm^ of 
cumulative liquid injected), which was difficult to do with a flowmeter or a differential 
detector, 

- to operate the measuring instrument with liquids of different viscosities without 
5 having to adapt the system for measuring the volume injected, which greatly improves 

the flexibility of the instrument. 

Figures 3A to 3C show examples of injected volume curves as a function of time 
for three different rocks. Significant variations of the filling curves are observed 
according to the nature of the rocks tested. The lower the permeability of the rocks, the 
10 slower the filling kinetics (more than 80 seconds for the least permeable rock shown, 
Fig.4A, and less than 2 seconds for the most permeable rock shown, Fig.4C). 

The goal of the next two stages is to obtain, fi-om the pressure or injected liquid 
volume measurements, an estimation of permeability K alone (if the porosity has been 
measured before as described above) or a combined estimation of permeability K and of 
15 porosity (4)). 

2) Modelling 

The drill cuttings are considered to be of homogeneous size and of spherical shape, 
and the gas is assumed to be perfect. The pressure drop due to the viscosity of the gas is 
disregarded in relation to that of the liquid, considering the difference between the 
. 20 viscosities. The residual gas trapped in the rock fi*agments after spontaneous imbibition 
of the liquid has the form of disconnected homogeneously distributed pockets. The 
capillary pressure is also considered to be negligible. 
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Considering the spherical shape of the fragments, we will consider a cap of 
thickness dr (Fig.2) and calculate the evolution of the pressure at the boundary of the 
rock particle when a flow of liquid q is injected. 

We consider that the total flow rate Q of the liquid injected is equitably divided 

5 among the N rock particles and that each one receives the flow rate q = %^ The gas law 

N 

allows to deduce the local gas saturation Sg from the moment that the pressure P : Sg = 
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p 

(Po is the pressure of the liquid) is known. In the cap, a material balance is 



performed on the liquid. The accumulation is equal to the difference between what goes 
in and what goes out. We deduce therefrom that : 



at 



Since So = (1 - 5^) = (1 - -B> we deduce that : 



P 



dSo _dSodP P^dP 



dt dP dt V ^°P^Jdt 
— K 

Besides, since = gradP^ (K is the permeability and fio the viscosity of the 

liquid), and since the capillary pressure can be considered to be negligible, so that we 
have Po = Pgss = P, the previous equation can be written in the following form : 

Mo P' dt 

15 It follows therefrom that : AP = . 

K P"" dt 
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We therefore obtain the conventional form of a diffusion type equation with^ 
however, a tenn factor of the accumulation due to the compressible nature of the 
gas. 

In spherical coordinates, the Laplacian is -j-— — I Finally, the equation to be 

r dr\ or J 

5 solved is written as follows : 



5 r 2 ap V dP ... 



With ^ = (2). 

K 



As it is injected, the liquid e?q)els the air in the free space between the drilling 
fragments and it enters the rock by spontaneous imbibition. Despite certain precautions, 
a certain volume of gas may remain outside because of the non-regular shape of the drill 
10 cuttings. This trapped volume (Vgp) plays a direct part in the general form of the 
pressure response and it has to be taken into account in the solution. 

A certain compressibility due to the experimental device also has to be taken into 
account. It comes from the chamber, from the lines as well as from the properties of the 
liquid. The equivalent compressibility observed is of the order of 0.0005 bar"^ (50 Pa). 

15 As the liquid used is saturated with gas at atmospheric pressure, dissolution 

phenomena occur when the pressure rises during measurement. These aspects are taken 
into account by introducing a diffusion parameter representative of the molecule 
exchanges at the level of the gas/liquid interfaces. 
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The diffusion equation is solved by means of the finite-difference method with an 
explicit scheme and by applying the boundary conditions in time P(r,0)=Patai and in 

space P(R,t)=Pext, and —(0,/) = 0 . 

or 

Within the context of the experimental procedure described above, the pressure 
5 evolution in the chamber containing the fi-agments, i.e. Pext, is directly measured. It is 
thus possible to directly solve the equation governing the evolution of the pressure in 
the rock explicitly. The evolution of the simulated volume of liquid injected as a 
function of time is then calculated by difference between the volume of gas initially 
present in the rock (at Patm) and the volume occupied by the gas at a time t during the 
10 test (during the compression stage). 

3) Adjustment of the model to the experimental results 

The model is implemented in a calculator such as computer 8 (see Fig. 1) in form of 
a software and it is inserted in an iterative optimization loop. The model is operated 
with a priori selected permeability K and residual gas saturation values, by imposing for 
15 porosity (<D) the value obtained experimentally, and the resulting simulated injected 
volume curve is compared with the experimental curve, by successive iterations, by 
changing the values of K and of the residual gas saturation in the model. The values that 
allow the best possible adjustment of the theoretical and experimental curves are found 
through the agency of a conventional Newtonian or gradient type optimization method. 

20 In cases where porosity (O) has not been previously measured experimentally, the 

model is operated with a priori selected permeability K, residual gas saturation and 
porosity values, and the resulting simulated injected volume curve is compared with the 
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experimental curve, by successive iterations, by changing the values of K, of porosity <1> 
and of the residual gas saturation in the model. Figure 1 1 shows that a satisfactory 
estimation of the porosity of the rocks tested is obtained. 

When the value of the previously measured porosity is known, only the 
5 permeability value K has to be adjusted in the modelling stage. The modelling operation 
is therefore faster. It is also more accurate, as can be clearly seen in Figure 12. 

Figures 6A to 6C show the comparison obtained by successive iterations between 
the theoretical and experimental curves for three rock types tested. In any case, a very 
good match is observed between the simulations and the experimental results. The 
10 results obtained by applying the method are quite comparable to those obtained in the 
laboratory after long conditioning periods using conventional methods for several rocks 
of different permeabilities as shown in Fig. 5. 

Programming this modelling within a code allows to calibrate the experiments by 
means of an optimization loop, which allows the corresponding value of K to be rapidly 
15 deduced. 



